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. Abstract 

Ifl ■ This paper proposes the density and characteristic functions of a general matrix 

(~| . quadratic form X*AX, when A = A*, X has a matrix multivariate elliptical dis- 

tribution and X* denotes the usual conjugate transpose of X. These results are 
obtained for real normed division algebras. With particular cases we obtained the 
density and characteristic functions of matrix quadratic forms for matrix multivari- 
ate normal, Pearson type VII, t and Cauchy distributions. 



■ 1 Introduction 

■ The distribution of a matrix quadratic form in multivariate normal sample has been studied 
VO I by diverse authors using zonal, Lagucrrc and Hayakawa polynomials w ith matrix a rgument, 
(-^ ■ in real and complex cases, see iKhatril (119661), iHavakawal (20071) and IShahl (|l970[ ). among 



others. It is important to observe that these results have been obtained for positive definite 
OA ■ matrix quadratic forms. 

In diverse cases, certain statistics are functions of a quadratic form or special types of 
it, and play a very important role in classical multivariate statistical analysis. 
_ By replacing the matrix multivariate normal distribution with a matrix multivariate 

■ elliptical Distribution, in classical multivariate analysis one obtains what is now termed 

I generalised multivariate analysis. By analogy, the distribution of a matrix quadratic form 

in a matrix multivariate elliptical sample plays an equally important role in generalised 
multivariate analysis, as it is now interesting to study the distribution of a quadratic form 
assuming a matrix multivariate elliptical distributio n. Some results in this conte xt have been 
obtained for particular matrix quadratic forms, see iFang and Anderson (Il990l Chapter II, 
pp. 137-200). 

In general, many results first described in statistical theory are then found in real case, 
and the version for complex case is subsequently studied. In terms of certain concepts and 
results derived from abstract algebra, it is possible to propose a unified means of addressing 
not only real and complex cases but also quaternion and octonion cases. 

In this paper we obtain the density and characteristic functions of a matrix quadratic 
form of a matrix multivariate elliptical distribution for real normed division algebras. Fur- 
thermore, these results are obtained when the matrix of the quadratic form is not necessarily 
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positive definite, which generalises most of the results presented in the literature in this con- 
text. This paper is structured as follows: Section [2] provides some definitions and notation 
on real normed division algebras, introducing the corresponding matrix multivariate ellip- 
tical distributions. Some results for Jacobians are proposed and two are obtained together 
with an extension of one of the basic properties of zonal polynomials. This is also valid for 
Jack polynomials for real normed division algebras, termed spherical functions for symmet- 
ric cones, and are also obtained. Section [3] then derives the matrix quadratic form density 
function, and as corollaries, some results for particular elliptical distributions are obtained. 
In Section 21 the characteristic function of a matrix quadratic form is obtained and some 
particular cases are studied. 

2 Preliminary results 

Let us introduce some notation and useful results. 



2.1 Notation and real normed division algebras 

A comprehensive discussion of real normed division algebras can be found in lBaezl(|2002l ). For 
convenience, we shall introduce some notations, although in general we adhere to standard 
notations. 

For our purposes, a vector space is always a finite-dimensional module over the field 
of real numbers. An algebra 5^ is a vector space that is equipped with a bilinear map 
m ■ S >^ -S ^ S termed multiplication and a nonzero element 1 G 5^ termed the unit such 
that m{l, a) — m{a, 1) = 1. As usual, we abbreviate m(a, b) = ah as ab. We do not assume 
^ to be associative. Given an algebra, we freely think of real numbers as elements of this 
algebra via the map cj i— >■ wl. 

An algebra is a division algebra if given a, 6 G 5^ with ab — 0, then either a = or 
6 = 0. Equivalently, 5^ is a division algebra if the operation of left and right multiplications 
by any nonzero element is invertible. A normed division algebra is an algebra 5 that is 
also a normed vector space with ||a6|| — ||a||||6||. This implies that 3^ is a division algebra 
and that ||1|| = 1. 

There are exactly four normed division a lgebra s : real numbers (3?), complex numbers 
(£), quaternions (io) and octonions (D), see iBaej (j2002l) . We take into account that 5R, 
£, f) and D are the only normed division algebras; moreover, they are the only alternative 
division algebras, and afidiyision algebras have a real dimension of 1,2,4 or 8, which is 
denoted by /?, see iBaej ( 2002 . Theorems 1, 2 and 3). In o ther branches of mathematics, the 
parameter a = 2/^ is used, see lEdelman and Raol (i2005h . 

Let „ be the linear space of all n x m matrices of rank m < n over ^ with m distinct 
positive singular values, where 5 denotes a real finite- dimensional normed division algebra. 
Let 5"^™ be the set of all n x m matrices over ^. The dimension of g^"^™ over 5R is /Smn. 

T 

Let A G 3^"^™, then A* = A denotes the usual conjugate transpose. 

The set of matrices Hi g J^x™ such that H*Hi = I™ is a manifold denoted „, termed 
the Stiefel manifold (Hi, also known as semi- orthogonal (/? = 1), se mi-unitary (0 = 2), 
semi-s ymylectic {13 — 4) and semi-exceptional type {j3 = 8) matrices, see Dray and Manogua 
()l999f) i. The dimension of „ over K is [f5mn — m{m — l)/3/2 — m]. In particular, V£ ^ 
with a dimension over 5R, [m(m -|- l)/3/2 — m], is the maximal compact subgroup il^{m) of 
£1^ ,^ and consists of all matrices H e 5^™''" such that H*H = I,„. Therefore, il'^(m) is the 
real orthogonal group 0{m) (/? = 1), the unitary group U{m) (/3 = 2), compact symplectic 
group Sp{m) (/3 = 4) or exceptional type matrices Oo{in) {j3 = 8), for = K, £, or £), 
respectively. 
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We denote by 6^ the real vector space of all S e such that S = S*. Let 

be the cone of positive definite matrices S S then is an open subset of 

Over 3fJ, consist of symmetric matrices; over £, Hermitian matrices; over Sj, 
quaternionic Hermitian matrices (also termed self-dual matrices) and over D, octonionic 
Hermitian matrices. Generically, the elements of shall be termed Hermitian matrices, 
irrespective of the nature of ^. The dimension of over is [m{m — l)/3 + 2]/2. 

Let be the diagonal subgroup of C^^^ consisting of all D £ g^mx™^ j) — diag(di, . . . , dm) 
and let 1u^^^ be the semi-upper-triangular subgroup of -C^ „ consisting of all T G ^"xm^ 
with > 0. 

For any matrix X G ^"x™^ denotes the matrix of differentials {dxij). Fi nally, we 
define the measure or volume element (dX) when X G 5-™^", S^,, 2)^, or V/^ „, see lPimitriul 
(120021 ). 

If X G g'"^™ then (dX) (the Lebesgue measure in g^"xm-j denotes the exterior product 
of the j3mn functionally independent variables 



(dX) — f\ f\ d-Xij where dxij — ^ dx[ 



4f. 

i=ij=i k=i 



(dS) 



If S G &t (or S G 2:^(m)) then {dS) (the Lebes gue measure in &(^ or in 1.^{mj) denotes 
the exterior product of the m{m + l)/3/2 functionally independent variables (or denotes the 
exterior product of the m(rn — l)/3/2 + n functionally independent variables, if su G 5ft for 
all i = 1, . . . , m) 

m /3 

i<j k^l 
m m 

t\dsu f\ t\dsf, ifs,, G5ft. 

i—1 i<j k—1 

The context generally establishes the conditions on the elements of S, that is, if Sij G 3?, 
^ G Si OT e D. It is considered that 

m (3 m rn /3 

{ds)^/\/\dsf^/\ds.,/\/\d^!;y 

i<j k—1 i—1 i<j k—1 

Observe, too, that for the Lebesgue measure (dS) defined thus, it is required that S G 
that is, S must be a non singular Hermitian matrix (Hermitian positive definite matrix). 

If A G S)^ then (rfA) (the Legesgue measure in Sf^) denotes the exterior product of the 
/3m functionally independent variables 

n p 

(dA) = /\/\dAf). 

i=l k=l 

If Hi G V£_„ then 

n m 

(H^dHi) = /\ /\ h*dh,. 

i=i j=i+i 

where H = (H1IH2) = (hi, . . . , hm|hm+i, . . ■ , h„) G ii'^{m). It can be proved that this 
differential form does not depend on the choice of the H2 matrix. Furthermore, (HirfHi) 
is invariant under the matrix transformations 

Hi^HiQ, Q Git'^(m) and Hi ^PHi, PGii'^(n), 
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and defines an invariant measure on the Stiefel manifold „, see Jamed ( 19541 ). 

When m = 1; vf „ defines the unit sphere in y. This is, of course, an {n — 1)/?- 
dimensional surface in 5^". When m = n and denotes Hi by H, (H*dH) is termed the Haar 
measure on 

The surface area or volume of the Stiefel manifold V£ „ is 



Vol(V/^,„ 



^2'^ j^inn P / 2 



(H^dHi) = ,, 



(1) 



and therefore 



(rfHi) 



Vol 



is the normalised invariant measure on V/^ „ and (dH), i.e., with (m = n), it defines the 
normalised Haar measure onil^(m). In ([1}, r^[a] denotes the multivariate Gamma function 
for the space 6^, and is defined by 

ri[a] = / etr{-A}|Ap-("-i)^/2-i(dA) 



IT 



m 

r„(m-l)/3/4-Qpj^_(^_^)^/2], 
i=l 



where etr(-) = exp(tr(-)) , | ■ | denotes the determinant and Re(a) > (m — l)/3/2, see 
Gross and Richardd (|l987t ). 



2.2 Jacobians 

In this section, we summarise diverse Jacobians with respect to Lebesgue measure i n term s 
of the parameter. For a detailed discuss ion of this and re lated issues, see iDimitriul (|2002l ). 
Edelman and Raol (|2005l ). lForresteij ([2009h and lXab^ (|l984l ). Two new Jacobians for singular 



linear transformations are proposed as well. 

Proposition 2.1. Let A £ £^ B G £^ „ and C g £^ „ be matrices of constants, Y and 
X G „ matrices of functionally independent variables such that Y = AXB + C. Then 

(dY) = |A*A|'3™/2|B*g|/3„/2(^x). (2) 

Proposition 2.2 (Singular value decomposition, SVD). Let X G £„„, such that X = 
ViDW* with Vi G Vf„^„, W G H'^(m) and D = diag(di, • • • , d™) G 2)^^', di > • • • > d,„ > 
. Then 



(dX) = 2-"V^ J] -Q^^; 

i—1 i<j 



'2 d^.fidT))iYldYi)iW*dW), (3) 



whe 



0, 




= 1; 


— m. 


/? 


= 2; 


-2m, 


/3 


= 4; 


—4m, 


/3 


= 8. 
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Proposition 2.3 ( Spectral decomposition). Let S 6 ^„i- Then the spectral decomposition 
can be written as S — WAW*, where W G il'^(m) and A = diag(Ai, . . . , Am) G 2)m; with 
Ai > • • • > Am > 0. Then 



(dS) 2-™7r^ ]J(A, - Aj )'^(dA)(W*dW), 



(4) 



where r is defined in Proposition 
Proposition 2.4. Lei X e £f„.„, and S = X*X e T/ien 

(dX) = 2-™|S|'3("-'"+i)/2-i((iS)(VtdVi), (5) 

Now, let £^ „(q) be the linear space of all n x m matrices of rank q < min(n, m) with q 
distinct singular values. In addition, observe that, if X G n{l)j can write X as 



Xi 



Xii Xi2 

qxq qxrn—q 

X21 Xll 

n — qxq n—qxm — q 



such that r(Xii) — q. This is equiv alent to the left and right pr oducts of matrix X with 
permutation matrices IIi and 112, see lGolub and Van LoanI (|l996l section 3.4.1, 1996), that 
is Xi = niXn2. Note that the exterior product of the elements from the differential matrix 
dX is not affected by the fact that we multiply X (right and left) by a permutati on matrix, 
that i s, (dXi) = (dfniXna) ) = (dX), sin ce Hi G ii'^(n) and U2 G il^(m), see iMuirhead 
(|l98l Section 2.1, 1982) and I JamesI (|l954h . Then, without loss of generality, (dX) shall be 



defined as the exterior product of the differentials dxy , such that Xij are mathematically 
independent. It is important to note that we shall have {nq + mq — q'^)(3 mathematically 
independent elements in the matrix X G C^^ nil)^ corresponding to the elements of Xn, X12 
and X21. Explicitly, 



(dX) = (dXn) A (dXi2) A (dX2i) = /\ /\ /\ dx\;> /\ /\ /\ dx 



(6) 



=1 j=i k=i 



= lj=q+l k=l 



Furthermore, 

(dX) = (PdXQ), for P G il^(n), Q G il^(m) (7) 

Observe that an explicit form for (dX) depends on the factorisation (base and coordinate 
set) employed to represent X or, that is, they depend on the measure factorisation of (dX). 
For example, by using the non-singular part of the decomposition in singular values and the 
non-singular part of the spectr al decomposit io n for X, then we can find a n expli c it form for 
(dX), which is not unique. seelKhatril (Il968h . lOfaz-Garcfa et al. I (|l997t). lUhlig I (|l994[) and 
Dfaz-Garcfa and Gutierrez- Jaime 1997t ). Alternativelv. an explicit form for (dX) can be 
found in terms of the QR decompositions. 

A singular random matrix X in „ (q) does not have a density with respect to Le besgue's 
measu r e in g"^™ but it d oes possess a density on a subspace M C g""^™; see iKhatri 
()l968l ). lRaol(ll973l p. 527'). lDfaz-Garcfa et al. I |l997l ). IUhhg I (|l994l) and fcramer I ()l999l p. 



297). Formally, X has a density with respect to Hausdorff's measure, w hich coincides with 
Lebesgu e's m easure, when the latter is defined o n the subspace Ai; s ee Billingsle^ (11986 



p. 2 47').IUhhg ( 19941 .iDfaz-Garcfa et ah I (119971 ). iDfaz-Garcfa and Gutierrez- Jaimed (|2005[ ) 



and Dfaz-Garcfa and Gutierrez- Jaimed (|2009 ) . 



Given X G „(g), and constant A G £^ p(c), and Y G ^(9) with c> q. We wish to 
determine the Jacobian of the transform Y = AX. Let us first consider the following case: 
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Theorem 2.1. LetX G £^j_„(7i), with A 6 £fi_p(r7.) constant andY e £,^p{n). IfY = AX, 
then 

n n 

{dY) = II a,(A)^'"(dX) = Yl A,(AA*)^"/2(dX) (8) 

1=1 i=l 

where Ai(M) anc? cri(M) are t/ie i-t/i non-null eigenvalue and singular value o/ M, respec- 
tively. 

Proof. Let A = HiDaQ* is the non-singular part of the SVD of A, where Hi G V^p, 
Da = diag(o-i(A), • • • ,cr„(A)), with (Tj(A) the z-th singular value of A and Q e il^(n). 
Furthermore, note that rank(Y) = rank (AX) = n. By differentiating Y = AX, we obtain 

dY = ArfX = HiDAQ*dX. 

Now, let H2 (a function of Hi) be such that H = (Hii H2) S il^(n), then 

WdY ^ r Sn HiDaQ^^X - f HtHiDAQ*rfX ^ _ f DaQ^^X 



H2 J \^ H*HiDAQ*rfX J \ 

as H^Hi = 0. From (0), (H*rfY) = (dY) and that (Q*dX) = (dX), then by Proposition 

n n 

{dY) = |DA|'^"(dX) = J|a,(A)^"(dX) = J| Ai(AA')''™/'(dX). □ 

1=1 i=l 

□ 

Remark 2.1. Note that, we can consider the QR decomposition instead of the SVD of 
matrix A in Theorem \TT\ That is A = HiT, where Hi G V^^p and T G 1u+^, see 



Diaz-Garcia and Gutierrez- Jaimez ( 19971 ). Alternatively to (O we have that 

n 

{dY) = l[ti"\dX). (9) 

1=1 

The proof is parallel to that given in Theorem 12. II Additionally note that, when n = p, ^ 
and © they agree. 

Theorem 2.2. Let X G £^,„((7), with A G £^_p(c) constant, and Y G £^,p((?), wit/i 
min(p, n) > c > g. IfY = AX, t/ien 

J|A,(ACC*A*)'^"/2 
(dY) = i^i^ (dX) (10) 

1=1 

w;/iere C G £^,„(g). 

Proo/. Let C G £^,„(g) such that X = CZ where Z G £(^.,(9) and let us denote R = AC. 
Then 

Y = AX = ACZ = RZ. 

Observing that rank (Y) = rank (RZ) = rank (Z) = q, from Theorem 12. II we have 

(dY) = [|A,(RR*)'^"/2(^2) = [|A,(ACC*A*)'5'"/2(dZ). (11) 

i=l 1=1 
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Now, X = BZ, again applying Theorem 12. 11 we obtain 

q 

i=l 

from which, substituting (dZ) = JlLi ^i(CC*)"'^"/2(^x) in (HH), we obtain the desired 
resuh. □ □ 

Remark 2.2. Without loss of generality, suppose that C g such that C*A*AC = 
diag(Ai(A*A), . . . , Aq(A*A)), then from Theorem [TUl 

(dY) = '[[cT.iAf'^idX) =f[X,{AA*f^/^{d-K) 

i=l i=l 



2.3 Elliptical distributions 

Now, the generalised matrix multivariate elli ptical distributions for real normed divi sion 

algebras it is introduced in current section, see iDi'az-G arcfa and Gutierrez- Jaimez ( 1997). A 

compr ehensive and systematic study can be found in[Fang an d ZhanJ i 19901 ) and Gupta, and Varga 
( 1993), for rea l case . Similarly, the elliptical vector distributions have been discussed by 



Micheas et al. for complex case. 



Definition 2.1. It is said that the random matrix Y £ „ has a matrix variate elliptical 
distribution, denoted as Y ~ £^^j^{fi,@,Tl,h), if its density with respect to the Lebesgue 
measure on y ^™ is 

C^{m,n) 



|X;|/3"/2|e|/3™/2 

where 



/i{/3tr[S-i(Y-M)*e-i(Y-/x)]}. (12) 




C'im,n) = "^fZr \ I . u^—'h{Pu-)du \ (13) 



and where 6 *pf,, S G «pf„ and ^ e „ are constant matrices. 

Observe that this class of matrix multivariate distributions includes normal, contami- 
nated normal, Pearson type II and VII, Kotz, Jensen-Logistic, power exponential and Bessel 
distributions, among others; these distributions have tails that are more or less weighted, 
and/or present a greater or smaller degree of kurtosis than the normal matrix multivariate 
distribution. In particular, observe that if in Definition [2TT] it is taken that h{u) = exp(— u/2), 
from (jl3p it can be readily seen that Cl^{m,n) = (27r^-i)-™"'^/2^ Hence, the density ob- 
tained is 



1 



(27r/3-l )™"'3/2 I S|'3«/2 I |/3m/2 



etr tr [Y^-\Y - m)*0-1(Y - /x)] | 



(14) 



which is termed the matrix multivariate normal distribution for real normed division algebras 
and is denoted as Y '-^ Af^^^{fi, 0, S). 

Similarly, observe that if in Definition 12.11 it is taken that h{u) = (1 + u/g) where 
s,g s,g > 0, s > /3mn/2; from (fT3| it can be seen that 



(7r<?/3-i)/3W2rf 



/3mn 
2 



7 



Therefore, the density is 



~2~ 



1 + ^ tr [S-i(Y - /x)*0-i(Y - /x) 

£;|/3n/2|0|/3m/2 ^ 9 



(15) 



which is termed the matrix multivariate Pearson type VII distribution for real normed 
division algebras. Observe that when s — {/3mn + g)/2 in (|15p. Y is said to have a matrix 
multivariate t distribution for real normed division algebras with g degrees of freedom. And 
in this case, if g = 1, then Y is said to have a matrix multivariate Cauchy distribution for 
real normed division algebras. 



2.4 Some results on integration 

First consider the following result proposed bv IConstantine and Muirhead ( 1976) for real 
case a nd extended to real normed division algebras by iDiaz-Garcia and Gutierrez- Jaimez 
( 1997t ). Observe that given a function /(H), H S il^(n) is possible to integrate over the last 
n — m (n > m) columns of H, the first m columns being fixed, and then to integrate over 
these m columns, that is 



Lemma 2.1. 



/(Hi,H2)(H*(iH) 

U'3(m) 



HiGV,, 



/(Hi, FM)(M*dM)(H*dHi), 

Meil'3(n-m) 



where H = 
those of Hi 



(HiiHa), Hi G 

, this is FF* = 



V,l„ and F 
I — HiHi. 



F(Hi) e V, 



n—m,n 



with columns orthogonal to 



Jack polynomials for real normed division algebras are also te r med s pherical functions 
of symmetric cones in the abstract algebra context, see ISawveii (|l997l ). In addition, in 
the statistical literature, they are termed real, complex, quaternion and octonion zona l 
polynomials , or, generically, general z onal polynomials, see I JamesI (|1964 ). iMuirheadl (|l982r ). 
Kabd (|l984l ). lRatnaraiah et all (|2005[ ) and lLi and Xud (|2009l ). 

Th e ori ginal ve r sion o f the follow result was stated for real and complex cas e s in JamesI 
( 1960) and James (1964'), respectively, and for quaternion case bv iLi and Xue ( 20091 ). see 
also iDiaz-Garcfa (,2009 ). Next these results are ex tended under mo re general conditions. 
This result was previously indirectly conjectured by iHavakawal (|2007i Remark on p. 194). 

Lemma 2.2. Let Xi e 6f„ and Xa G 6f„ and H e ^'^(to). Then 



L 



HGU'5(n 



(XiHX,H*)(dH) = '^'^^^^^''^^ 



(16) 



where C^{-), denotes the Jack polynomials fsee lSawvei 1 199't ) and \Diaz- Garcia 1 200!^ )). r 
— rank of X2 and (dH) is the normalised invariant measure on ii^ (n) . 

Proof. As well as taking into account that this ne w version pr oposed holds fo r real normed 



m pr e . 

division algebras, in the original result proposed James (|l960i) and JamesI fl964), X2 is 
assumed to be positive definite, and in ou r case X2 i s assu med positive semidefinite. Then, 
the proof is the same as that given by I Muirhead ( 1982 . Theorem 7.2.5, pp. 243-244), 
simply observin g that the differential operator Ay„ is also invariant when X2 is positive 
semidefinite, see iDiaz-Garcia and Caro-Lopera (1200^ . Hence, we take 



X, = 



Ir 







r xr—rn 



r — m X rri r—m xr — ni 
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where r is the rank of X2, the conclusion is obtained observing that 

a(A) = C«(^^ (17) 

iKushneij (|l985l . Equation (3.6), p. 88), see also iMuirheadI (|l982l Corollary 7.2.4(i), p. 
236). □ 

□ 

Theorem 2.3. Let Xi G anrf X2 £ 6f^, and H = (HiiHa) G il'3(n), wit/i Hi £ V™,„. 
T/ien 

^ , Cf(Xi)C,^(X2) 

cf(I.) ^ ^ 

where r — rank of X2 . 

Proof. First, denote by J the integral in the left side of ([TSl) . Then, for K e U'^(n — m) 
J = / / Cf(XiHiX2Ht)(dK)(HtdHi) 

where (dK) denotes the Haar probability measure on il^ [n — m). This is by ([T]) 
(^^) = TTli^ ^(K*rfK), from where, / (dK) = 1. 

V0l(ilP(n - m)) JKGU/3(n-m) 



/ (XiHiX2HJ)(HtdHi) - Vol(V, 



Then, for H = [HiiHa] G il'3(n) and by Lemma O 

/ Cf(XiHiX2Ht)(H*dH) 



Vol(il'3(n-m)) yneii^H 
Now, normalising the measure (H*(iH) and noting that 

Vol(il^(n)) 
Vol(ii'3(n-m)) 

we obtain: 



= Vol(Vi,„) 



J = Vol(Vf„,j/ C,^(XiHiX2HI)(dH) 

"'Hell'^(n) 



Finally, observe that 



Cf(XiHiX2HD = Cf (xi(Hi;H2) ( "^^r" US 



I V H 

ii—nixin n — mxn — m 



the desired result is obtained from (fTTl) and Lemma [2?2l □ □ 

A basic integral property is cited below. For this purpose, we utilise the complexification 
6^a^ = + iet of 6f„. That is, S^;'^ consists of aU matrices X G (5C)™x™ of the form 
Z = X + iY, with X, Y G 6f^,. We refer to X = Re(Z) and Y = Im(Z) as the real and 
imaginary parts of Z, respectively. The generalised right half-plane $ = + i&^ in 
consists of ah Z G 6^*^ such that Re(Z) G see iDiaz-Garcia (,2009, ). 
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Theorem 2.4. Let Z E ^ and U e 6f„. Assume 7 = J^^^ii f{z)z'''^-''-'^dz < 00. Then 
[ /(trXZ)|X|'^-("-i)^/2-i^^ (XU) (dX) 

i [am + k\ 

where 1} = J^^^^ f{z)z°-™'^''~^dz, Re(a) > (m — l)/3/2 — fcm a'^c^ = (fei, ■ • ■ , fcm) k = 

fci + • • • + fcm, aee \Diaz- Garctd l200^ ). And where [a]^ denotes the generalised Pochhammer 
symbol of weight k, defined as 

m 

{a]t = ilia - (z - l)/3/2),^ ^ , 

where Re(a) > (m — l)/3/2 — km and 

(a)j = a{a + 1) • ■ • (a + i - 1), 
is the standard Pochhammer symbol. 

Finally, consider the extension of the Wishart's integral for real normed division algebras. 
Theorem 2.5. Let Y e £f„.„. Then 

/ /(Y*Y)d(Y) = |R|^("-"'+^)/^-V(R). 

Proof. By the Proposition [2]4] 

(dY) = 2-"|R|'^("-"+i)/2-i(rfR)(VtdVi), 

with Vi e V,*^ „. Hence 

/ /(Y*Y)(i(Y) = 2-"/" /" |R|'^("-'"+i)/2-V(R)(dR)(V*dVi). 
Jy*y=r -/Re'^^m "'Viev^,„ 

from which the conclusion is reached. □ □ 

3 Density function 

In this section we find the density function of a general matrix quadratic form of a matrix 
multivariate elliptical distribution for real normed division algebras. 

Theorem 3.1. Assume that X - x-m(0, 0, S, h) and define W = X*AX G ©f„ of rank 
r, with A G (3^ 0/ ranA; r < m < n. Then the density function of W for a real normed 
division algebra, is 

[/3n/2] I S |/5"/2 |0 |/3m/2 I A[/3™/2 ^ ~fc! C^(Ir) 

where A"*" is t/ie Moore-Penrose inverse of A; A = PiAP* is the non-singular part of the 
spectral decomposition of A, whit Pi G Vrm and A = diag(Ai (A), . . . , Ar(A)), Ai(A) > 
• ■ • > Ar(A) > 0; and h^''\-) is the kth derivative of h. 
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Proof. The density function of W is 



[ h[l3tT S-iX*©-iX] (dX) 

JX*AX=W 



|S|/3«/2|0|^W2 7x.AX=W 

Let us now consider the transformation Y = A^/^X, such that (A^/^)^ = A, then, X 
A+^/^Y, where A+^/^ is the Moore-Penrose inverse of A^/^. Since, by Remark [221 



(dX) = =[]A,(A+i/2A+i/2*)/3W2(rfY) = []A,(A+)^"/2(^Y) 

r 

= ]^A,(A)-^"/2(^Y) = |A|-'3"/2(dY) 

where r is the rank of A; A = Pi AP* is the nonsingular part of the spectral decomposition 
of A, whit Pi G V^,„ and A = diag(Ai(A), . . . , Ar(A)), Ai(A) > • • • > Xr{A) > 0. Hence, 



Cl 



Y'Y=W 



where 



Cl = 



/3 tr Y* A+i/^e^i A+i/2 Y 



(dY), 



|2|/3n/2|0|/3m/2|A|/3m/2- 

Since A+^/^0^^ A+^/^ e 6^, the integral is invariant under the matrix transformation 

A+1/20-1A+1/2 ^ h*A+i/20-1A+i/2h, H e il^(n), 
and the integration with respect to H on il^ (n) . Therefore, 

/3trS-iY*H*A+i/20-iA+i/2HY 



Cl 



(dH)(dY), (21) 



where (dH) is the normalised invariant Haar measure. Let us now assume that h can be 
expanded in series of power, that is 

^ kl 

k=0 

Hence, recalling that 

5^Cf(X) = (tr(X))^ 

K 

Diaz-Garcia ( 2009( ). the /Hgu3(„) ^['K*^^) expression in (PT|) is 



see 



yhi^ r f/?trS-iY*H*A+i/20-iA+i/2HY')'(dH), 



k=0 

oc 



y y ^^—^ [ a f Y*H* A+i/20-1 A+1/2HY') (dH), 



and from Theorem 



/tW(O) q (/3YS-1Y*) (A+V20-1A+1/2) 
fc! 



fc = K 



/iW(0) C« (/3S-iY*Y) q (Q-1A+) 
fc! 



fe = K 
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Then, substituting in (PT|) . we have that the density function of W is 

Finally, the desired result is obtained from Theorem 12.51 □ □ 

Corollary 3.1. Assume that X has a matrix multivariate normal distribution for real 
normed division algebras. Let W = X*AX, with A G 6^ of rank r < m < n, A = 
PiAP* is the non-singular part of the spectral decomposition of A, whit Pi G V,*^^ and 
A = diag(Ai(A), . . . , Ar(A)), Ai(A) > ••■ > Xr{A) > 0. Then the density of W for real 
normed division algebras is 

Proof. This follows from (|22p noting that for the normal case, h{u) = exp{— m/2} and 
C^(to,7i) = (27r//3)-'5"™/2. □ □ 

In real case, this result has been f ound by iKhatri (1966) and iHavakawal ( 20071 ) in 
terms of zon al polynomials, by IShahl (Il970l) in terms of Laguerre polynomials and by 



Gupta, and Nagar ( 2000.) in terms of Ha yakawa polyno mials. All these results were proposed 
when A is definite positive. In addition. iKhatril (1966) obtain these results for complex case 
too. 

Corollary 3.2. Suppose that X has a matrix multivariate Pearson type VII distribution 
for real normed division algebras. Let W = X*AX, with A € (3^ of rank r < m < n, 
A = PiAP* is the nonsingular part of the spectral decomposition of A, whit Pi G and 
A = diag(Ai(A), . . . , Ar(A)), Ai(A) > ••• > At.(A) > 0. Then the density of W for real 
normed division algebras is 

r^[/3n/2]|S|'3"/2|0|/3m/2|A|/3™/2 ^ fc! (7^(1^) 

where 



(7rg;9-i)/3™"/2rf [s - Pmn/2] ' 

(s)j, s{s + 1) • • ■ (s + fc — 1), is the standard Pochhammer symbol, and s, g G 5R, s, g > 0, 
and 8 > mn/2. 

Where s, g £ ^, s, g > 0, s > Pmn/2. 

Proof. In this case we have 

C'^{m,n) = , and hlu) ^ (l + ul qV', 

(7rg/3~i)^"W2rf [s-/3mn/2] 

s,g £ 3?, s,g > 0, s > /3mn/2, see lDiaz-Garcia and Gutierrez- Jaime j (|l997l) . Then 

From where the desired result is follows: □ □ 
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4 Characteristic function 

In this section we derived the characteristic function of a general matrix quadratic form of 
a matrix muhivariate elhptical distribution for real normed division algebras. 

Theorem 4.1. Assume that X ~ ^^(O, 0, S, h) and define W = X*AX G (3^ of rank 
r, with A G o/ ranfc r < m < n. Then the characteristic function o/ W for a real normed 
division algebra, is 

« rf[/3TOn/2 + fc] fc! CM ^ ' 

where [a]f? is defined in Theorem \2.4\ and 

Proof. We have that 



^w(S) = Ew(ctr{*WS}) 

= Ex(ctr{iX*AXS}) 



/ etr{iX*AXS}/x(X)(dX) 



Considering the transformation Y = ^/^XS -'^Z^, where (B^/^)^ = B, then by Propo- 
sition EIH (dX) = |0|"/2|S|"/2(£;y). Hence from ([12]) the characteristic function of W 
is 

C^{m,n) [ etr{iSi/2Y*0i/2A0i/2YSi/2s}/i(^trY*Y)(dY) (24) 

Now, let Y = ViDW* the singular value decomposition of Y. Then, by Proposition [521 

m m 

(dY) = 2-'"7r^]Jdf<""™+^^"^[|(d2 -d|)'^(dD)(VtdVi)(W*dW). 

i—l i<j 

Therefore, the integral in is 

2-'"7r^ /" II etr{iSi/2WDV*0i/2A0i/2ViDW*Si/2S} 

m m 

X /i(/3trWD2w*) J|df^"""+'^"'[|(d2 _d2)'3(dD)(VJ(iVi)(W*(iW) (25) 

4=1 i<i 

Now, note that the integral on V,fj „ in ([55)1 that is 

/" etr{iSi/2wDVt0i/2A0i/2ViDW*Si/2S}(VJdVi) 

"'Viev,1,„ 

can be expressed as 

11Y.T\ , C^f (^Si/2WDVt0i/2A0i/2ViDW*Si/2s) (V^dVi) 

fc=i K ■ "'Viev^,„ 
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which by Theorem 12.31 is 

k=l K, K.Ly^yi-r) 

Hence, substituting in (pS)) we obtain 



^ ^ A;!Ck (Ir) "'weu/3(m) JDes 

m m 

X /i(/3trWD2w*)[]df'"-"+')-'n(^^' -4)^(^^^)(W*dW). (26) 

i— 1 z<j 

Now, considering the transformation V = WD^W*, then by Proposition 12.31 noting that 
A = D2, and (dA) = 2™|D|(dD), 

m m 

(dV) = W{<f^ - d2)/3 J]d,(dD)(W*dW), 

i<.j 2—1 

Thus, by (fM|) and the characteristic function of W is 



C^(m,n)2-™Vol(V^, 



fe=i K fc!C^(Ir) -'ve'p!?, ^ 



X /i(/3trV)|V|'^(""™+i)/2-i(dV). (27) 
FinaUy the desired resuh is obtained integrating ([27| using Theorem [231 D D 



Corollary 4.1. Assume that X /las a matrix multivariate normal distribution for real 
normed division algebras, and define W = X*AX S 6(^j o/ ranfc with A G o/ rank 
r < m < n. Then the characteristic function o/ W for a real normed division algebra, is 

! Cfs- ( I?^ ) 

Proof. Observe that, for the normal case, h[u) = exp{—u/2} and C'^(m, n) — (27r//3)~''™"/^. 
Then in this case 



?9 = / ^ /(z)z'''""/2+fc-i^^ 

cxp{-z/2}z''""/2+/c-i^^ ^ 2^™"/2+fc-i+/Jp/^[^^^/2 + A:]. 



Then the result follows. □ □ 

In the real ca se, i.e . /? = 1, this result was obtained bv lKhatril ()l966l) in terms of zonal 
polynomials y bv lShah ( 1970l ) in terms of Laguerre polynomials with matrix argument when 
A is a positive definite matrix. 

Corollary 4.2. Assume that X has a matrix multivariate t distribution for real normed 
division algebras, and define W — X* AX G of rank r, with A G of rank r < m < n. 
Then the characteristic function of W for a real normed division algebra, is 

V V C^(QA)C.(^g/?SS) 

^^4"fc! (s-/3mn/2-A:)fc C«(I.) ' ^ ^ 

where s,g G s,(? > 0, s > /3mn/2. 
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Proof. In this case we have 

C^(m,n) = , and h(u)^(l + u/gy 

(7r,g;3-i)/3"W2rf [s - /3mn/2] 

with s, 5 G Ji, s, 5 > 0, s > f3mn/2. Then 



§ = I /(z)z'''""/2+fe-i^^ 

{1 + z/gy z^"''"^^+''-^dz 



[/3mn/2 + A:]rf [s - l3mn/2 - k] 



_ fimn/2+k-l+l3 

Finally, recalling that rf [a] = [a + fc]/(a)fe for a + fc > 0, k an integer number. Then, 
taking a = s — Pmn /2 — fc 

[s - l3mn/2 - k] = ^ ^ ' ' ^ 



(s — l5mn/2 — fc)^ ' 

the desired result is obtained. □ □ 

In corollaries 13.21 and 14.21 observe that, when s = {f3mn + g)/2, W is said to have a 
matrix quadratic form t distribution for real normed division algebras with g degrees of 
freedom. And in this case, if = 1, then W is said to have a matrix quadratic form Cauchy 
distribution for real normed division algebras. 



Conclusions 

Any reader interested in a particular case - real, complex, quaternions or octonions - need 
simply take the part icular value of /3 in order to obtain the results desired. Furthermore, 
as is asseverated by iKabd (11984 . our results can be extended to hypercomplex cases, by 



simply substituting /3 by 2/3, obtaining the complex, bicomplex, biquaternion and biocto- 
nion (or sedenionic) cases. Observe that, alternatively to use the concepts and notation 
of real normed division algebras, it is possible to use the concepts and notation of simple 
Jordan algebras noting that one particular algebra more is contained. As was mentioned, 
there are four real norm ed division algebras and five Euclidean simple Jordan algebras, see 
Casalis and Letac (Il996l) . 



In summary, the density and characteristics functions of a matrix quadratics forms of 
matrix multivariate elliptical distribution are found under a unified approach that allows the 
simultaneous study of the real, complex, quaternion and octonion cases, generically termed 
distributions for real normed division algebras. In particular these results were particularised 
for matrix quadratics forms of a matrix multivariate normal, Pearson type VII, t and Cauchy 
distributions for real normed division algebras. 

Finally, note that with particular cases of the results obtained in this paper, one can 
obtain many of the results published in the literature. Thus for example, if in Corollarv l4.11 
A is a idempotent matrix and = 1, we obtained 

^^{h) - 2 p 2^2^— k\ CM ■ 

= 2''-i/3''""/2|I-2z/3SS)r^"/2^ 
and for real case V'w(S) = |I - 2iSS)|-"/2_ 
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